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There is compelling astrophysical evidence for the ex- 
istence of dark matter [l|, |2|, which amounts to about 
one quarter of the total energy density in the Universe. 
Yet, there is no experimental information on the parti- 
cle composition of dark matter [2|, |3j. A class of new 
physics models j4|, motivated by astro-particle observa- 
tions [5J, Iff , predicts a light component of the dark mat- 
ter spectrum. The bottomonium system of T states is an 
ideal environment to explore these models. Transitions 
r(3S) -» TT+iT-r{lS) and Y(2S) -> Tr+Tr-r{lS) offer a 
way to cleanly detect the production of T(1S) mesons, 
and enable searches for invisible or nearly invisible decays 



of the T(1S) 7}. Such decays would be a tell-tale sign of 
low-mass, weakly-interacting dark matter particles. 

The Standard Model process T(1S) — > ^vv is not ob- 
servable at the present experimental sensitivity |8j. An 
observation of T decays with significant missing energy 
would be a sign of new physics, and could shed light on 
the spectrum of dark matter particles x- The branch- 
ing fraction (BF) B(T(1S) — > xx) 1S estimated to be 
as large as (4-18)xl0- 4 [10, while B(T(1S) ->• jxx) 
is suppressed by 0(a), and the range 10~ 5 -10 -4 is ex- 
pected [8jj. 

The decays T(1S) — >• 7 + invisible might also proceed 



through Wilczek production [10| of an on-shell scalar 
state A : T(1S) -> 7A , A -> invisible. Such low-mass 
Higgs states appear in several extensions of the Standard 
Model (ll| . Constraining the low-mass Higgs sector is 
important for understanding the Higgs discovery reach 
of high-energy colliders [12]. The BF for Y(1S) -► 7A 
is predicted to be as large as 5 x 10~ 4 , depending on 
m^o and couplings [13]. If there is also a low- mass neu- 
tralino with mass m x < itiaoJ2, the decays of A would 
be predominantly invisible [l4j j . 

For multibody T(1S) —> jxx decays, the current 90% 
confidence level (C.L.) BF upper limit, based on a data 
sample of- 10 6 T(15) decays, is of order 10" 3 0. The 
limit on two-body T(1S) — >• 7 + V, X — > invisible decays 
is B{T{1S) -+7 + I)<3x 10- 5 for m x < 7.2 GeV [|. 
The limit on invisible decays of T(1S) is B(T(1S) — > 
XX) < 3.0 x 10- 4 0- 

This Letter describes a high-statistics, low-background 
search for decays T(1S) — > 7 + invisible, characterized by 
a single energetic photon and a large amount of miss- 
ing energy and momentum. This is the first search of 
this kind to use the Y(1S) mesons produced in dipion 
T(2S) -> ir + Tr~Y(lS) transitions. We search for both 
resonant two-body decays T(1S) —> -fA°, A° — > invisible, 
and nonresonant three-body processes T(1S) — > 7XX- 
For the resonant process, we assume that the decay width 
of the A resonance is negligible compared to the exper- 
imental resolution [16J . We further assume that both 
the ^4° and x particles have zero spin. The decays 
T(1S) — > 7XX are modeled with phase-space energy and 
angular distributions, which corresponds to S-wave cou- 
pling between the bb and XX- 

The analysis is based on a sample corresponding to 
an integrated luminosity of 14.4 fb~ collected on the 
T(2S) resonance with the BABAR detector at the PEP-II 
asymmetric-energy e + e~ collider at the SLAC National 
Accelerator Laboratory. This sample corresponds to 
(98.3±0.9) x 10 6 T(2S) decays. We also employ a sample 
of 28 fb~ accumulated on the Y(3S) resonance (Y(3S) 
sample) for studies of the continuum backgrounds. Both 
Y(3S) -> ir+ir-Y(2S) and Y(3S) -> 7r + 7i-r(lS) decays 
produce a dipion system that is kinematically distinct 
from the Y(2S) — > tt + tt~Y(1S) transition. Hence, the 
T(3S) events passing our selection form a pure high- 
statistics continuum QED sample. For selection opti- 
mization, we also use 1.4 fb~ and 2.4 fb~ datasets col- 
lected about 30 MeV below the Y (25) and Y(3S) reso- 
nances, respectively (off-peak samples). The BABAR de- 
tector, including the tracking and particle identification 
systems, the electromagnetic calorimeter (EMC), and the 
Instrumented Flux Return (IFR), is described in detail 
elsewhere (T7I [l8J] . 

Detection of low-multiplicity events requires dedicated 
trigger and filter lines. First, the hardware-based Level- 1 
(LI) trigger accepts single-photon events if they contain 
at least one EMC cluster with energy above 800 MeV. A 



collection of LI trigger patterns based on drift chamber 
information selects a pair of low-momentum pions. Sec- 
ond, a software-based Level-3 (L3) trigger accepts events 
with a single EMC cluster with the center-of-mass (CM) 
energy E* > 1 GeV [19j, if there is no charged track 
with transverse momentum pr > 0.25 GeV originating 
from the e + e~ interaction region. Complementary to 
this, a track-based L3 trigger accepts events that have 
at least one track with pr > 0.2 GeV. Third, an of- 
fline filter accepts events that have exactly one photon 
with energy E* > 1 GeV, and no tracks with momentum 
p* > 0.5 GeV. A nearly independent filter accepts events 
with two tracks of opposite charge, which form a dipion 
candidate with recoil mass (defined below) between 9.35 
and 9.60 GeV. 

The analysis in the low-mass region m^o < 8 GeV 
(m x < 4 GeV), which corresponds to photon energies 
E* > 1.1 GeV, requires the single-photon or the dipion 
trigger/filter selection to be satisfied; the trigger/filter 
efficiency for signal is 83%. In the high-mass region, 
7.5 < m A o < 9.2 GeV (3.5 < m x < 4.5 GeV), we only 
accept events selected with the dipion trigger/filter, since 
a significant fraction of this region lies below the energy 
threshold for the single-photon selection. This selection 
has an efficiency of 12.5% for signal events. 

We select events with exactly two oppositely-charged 
tracks and a single energetic photon with E* > 0.15 GeV 
in the central part of the EMC (-0.73 < costf* < 0.68). 
Additional photons with E* < 0.12 GeV can be present 
so long as their summed laboratory energy is less than 
0.14 GeV. We require that both pions be positively 
identified with 85-98% efficiency for real pions, and a 
misidentification rate of < 5% for low-momentum elec- 
trons and < 1% for kaons and protons. The pion can- 
didates are required to form a vertex with x^tx < 20 (1 
degree of freedom) displaced in the transverse plane by 
at most 2 mm from the e + e~ interaction region. The 
transverse momentum of the pion pair is required to sat- 
isfy pttttt < 0.5 GeV, and we reject events if any track 
has p* > 1 GeV. 

We further reduce the background by combining sev- 
eral kinematic variables of the dipion system [7| into 
a multilayer perceptron neural network discriminant 
(NN) [20j . The NN is trained with a sample of simulated 
signal events Y(1S) —> 7XX (m x = 0) and an off-peak 
sample for background; the NN assigns a value N close 
to +1 for signal and close to —1 for background. We 
require M > 0.65 in the low-mass region. This selection 
has an efficiency of 87% for signal and rejects 96% of the 
continuum background. In the high-mass region we re- 
quire N > 0.89 (73% signal efficiency, 98% continuum 
rejection). 

Two additional requirements are applied to reduce spe- 
cific background contributions. Neutral hadrons from the 
radiative decays Y(1S) -> "fK° L K° L and Y(1S) -> "/rin 
may not be detected in the EMC. We remove 90% of 



these background events by requiring that there be no 
IFR cluster within a range of 20° of azimuthal angle {<f>) 
opposite the primary photon (IFR veto). This selection 
is applied for m^a < 4 GeV and m x < 2 GeV, since the 
hadronic final states in radiative Y(1S) decays are ob- 
served to have low invariant mass [21] . 

For the high-mass range we suppress contamination 
from electron bremsstrahlung by rejecting events if the 
photon and one of the tracks are closer than 14° in <\>. In 
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addition, the two-photon process e + e 
e + e~7/', rj' — > "fir + Tr~, in which the e + e~ pair escapes 
detection along the beam axis and the two pions satisfy 
our selection criteria, produces photons in a narrow en- 
ergy range 0.25 < E* < 0.45 GeV. We take advantage of 
the small transverse momentum of the rj' and reject over 
half of these events by requiring the primary photon and 
dipion system to be separated by at most A0 = 160°. 
The signal efficiency for this requirement is 88%. 

The selection criteria are chosen to maximize e/(1.5 + 
\/B) 22]. where e is the selection efficiency for m x = 
and B is the expected background yield. The signal ef- 
ficiency varies between 2 and 11%, and is lowest at the 
highest masses (lowest photon energy). The backgrounds 
can be classified into three categories: continuum back- 



grounds from QED processes e + e" 



— > 77T+7T 



with 



particles escaping detection, radiative leptonic decays 
T(IS) — > j£ + £~, where leptons £ = e,/i,T are not de- 
tected, and peaking backgrounds from radiative hadronic 
decays and two-photon rf production. 

We extract the yield of signal events as a function of 
171 a" ( m x) m ^ ne interval < rn^o < 9.2 GeV (0 < m x < 
4.5 GeV) by performing a series of unbinned extended 
maximum likelihood scans in steps of to^o (m x ). We use 
two kinematic variables: the dipion recoil mass M rcco ;i 
and the missing mass squared M 
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(1) 
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where E* v is the CM energy of the dipion system, and V 
is the four-momentum. The two-dimensional likelihood 
function is computed for observables (M rccoi \, MJ.) over 
the range 9.44 < M rccoil < 9.48 GeV and -10 < M x < 
68 GeV 2 (low-mass region) and 40 < M| < 84.5 GeV 2 
(high-mass region). It contains contributions from signal, 
continuum background, radiative leptonic T(IS) back- 
ground, and peaking backgrounds, as described below. 
We search for the A in mass steps equivalent to half 
the mass resolution a(rriA°)- We sample a total of 196 
points in the low-mass < m^o < 8 GeV range, and 
146 points in the high-mass range 7.5 < m^o < 9.2 GeV. 
For the T(1S) — > "fxX search, we use 17 values of m x 
over < m x < 4.5 GeV. For each m^o (m x ) value, we 
compute the value of the negative log-likelihood NLL = 
— ln£(iV s i g ) in steps of the signal yield iV s ig > while 
minimizing NLL with respect to the background yields 



-Vcont (continuum), -/Vi ept (T(1S) — > -f£ + £~), and, where 
appropriate, iVhadr (radiative hadronic background) or 
Njji (two- photon rj' background) . If the minimum of NLL 
occurs for iV s i g > 0, we compute the raw statistical sig- 
nificance of a particular fit as S — y / 2\og(£/£ ), where 
Co is the value of the likelihood for iV s i g = 0. For small 
5, we integrate C(N s i g ) with uniform prior over 7V s i g > 
to compute the 90% C.L. Bayesian upper limits. In the 
range 7.5 < m^o < 8 GeV and 3.5 < m x < 4 GeV where 
the low-mass and high-mass selections overlap, we add 
NLLs from both datasets, ignoring a small (3%) correla- 
tion. This likelihood scan procedure is designed to handle 
samples with a very small number of events in the signal 
region. 



We use signal Monte Carlo (MC) samples [23J, |24| 
Y(1S) -> jA° and Y(IS) -> -yxx generated at 17 val- 
ues of rriAo over a broad range < m^a < 9.2 GeV and 
at 17 values of m x over < m x < 4.5 GeV to deter- 
mine the signal distributions in M\ and selection effi- 
ciencies. We then interpolate these distributions and ef- 
ficiencies. The signal probability density function (PDF) 
in M\ is described by a Crystal Ball (CB) function [25[ 
(T(1S) —> jA°) or a resolution-smeared phase-space 
function (T(1S) — > 7xx)- The resolution in M\ is 
dominated by the photon energy resolution, and varies 
monotonically from 1 GeV 2 at low m^o to 0.2 GeV 2 at 
m A a = 9.2 GeV. We correct the signal PDF in M x 
for the difference between the photon energy resolution 
function in data and simulation using a high-statistics 
e + e~ — > 77 sample. We determine the signal distribu- 
tion in M reco ;i, as well as that of background containing 
real T(1S) decays, from a large data sample of events 
T(IS) —t /J. + n~ ■ This PDF is modeled as a sum of two 
CB functions with common mean, a common resolution 
cr(M reco ji) Ki 2MeV, and two opposite-side tails. 

We describe the M| PDF of the radiative T{1S) -^ 
"f£ + £~ background by an exponential function, and de- 
termine the exponent from a fit to the distribution of 
M x in a T(1S) — > j£ + £~ data sample in which the two 
stable leptons (e or n) are fully reconstructed. Before 
the fit, this sample is re- weighted by the probability as a 
function of M x that neither lepton is observed. 

The continuum M x PDF is described by a function 
that has a resolution-smeared phase-space component at 
low M x , and an exponential rise at high M x . For the 
low-mass selection (—10 < M x < 68 GeV 2 ), we deter- 
mine this PDF from a fit to the T(3S) data sample. For 
the high-mass region (40 < M x < 84.5 GeV 2 ), we deter- 
mine this PDF, as well as the M x PDF of the peaking 
77' background, from a fit to the T(2S) data sample se- 
lected with the NN requirement TV < 0. The M rcco ii PDF 
is determined from a fit to the T(3S) data sample. 

The contribution from the radiative hadronic back- 
grounds is estimated from the measurement of T(1S) — > 
jh + h~ spectra [21] . We assume isospin symmetry to re- 
late B(T(1S) -> -yK+K-) to B{T(IS) -> jK° L K° L ), and 



B(T(1S) -> 7 pp) to B(T(15) -)■ Tnn). A small addi- 
tional contribution arises from T(IS') — > 77r + 7r~ events 
in which the pions escape detection. We expect A^adr = 
6.6 ±1.1 radiative hadronic events (without IFR veto), 
dominated by Y(IS) -> "/K° L K° L , or N^° = 1.02 ± 0.14 
events (with IFR veto). We describe the M\ distribution 
of these events with a combination of CB functions, using 
the measured spectrum of T(1S) — > "fh + h~ events [21 1 . 

The largest systematic uncertainty is on the recon- 
struction efficiency, which includes the trigger/filter ef- 
ficiency (etrig), and photon (e 7 ) and dipion (ew) recon- 
struction and selection efficiencies. We measure the prod- 
uct ew x Nrns), where N^as) is the number of pro- 
duced Y(1S) mesons, with a clean high-statistics sample 
of the Y(1S) —> ^ + M~ decays. The uncertainty (2.1%) 
is dominated by B(T(1S) —} /i + M~) (2%) [3J] and a small 
selection uncertainty for the /i + /i~ final state. We mea- 
sure e-y in an e + e~ — > 77 sample in which one of the 
photons converts into an e + e~ pair in the detector ma- 
terial (1.8% uncertainty). The trigger efficiency e t ri g is 
measured in unbiased random samples of events that by- 
pass the trigger/filter selection. This uncertainty is small 
for the single-photon triggers (0.4%), but is statistically 
limited for the dipion triggers (8%). In the low- mass re- 
gion, we take into account the anti-correlation between 
single-photon and dipion trigger efficiencies in L3; the 
uncertainty for the combination of the triggers is 1.2%. 

We account for additional uncertainties associated 
with the signal and background PDFs, and the predicted 
number of radiative hadronic events iVh a< j r , including 
PDF parameter correlations. These uncertainties do not 
scale with the signal yield, but are found to be small. 
We also test for possible biases in the fitted value of the 
signal yield with a large ensemble of pseudo-experiments. 
The biases are consistent with zero for all values of m^o 
and m x , and we assign an uncertainty of 0.25 events. 

As a first step in the likelihood scan, we perform fits 

0. 



to the low-mass and high-mass regions with N s { 
The free parameters in the fit are iVcont 1 Mept j an d -/Vhadr 
(low-mass region), and N cont , N\ ept , and N n > (high-mass 
region). The results of the fits are shown in Fig. [TJ We 
observe no significant deviations from the background- 
only hypothesis. We find JV hadr = 8.7^3;° ± 0.8 (without 
IFR veto) with a significance of 3.5er, including system- 
atic uncertainties. 

We then proceed to perform the likelihood scans as a 
function of N s { g in steps of m^o and m x . In the scan, the 
contribution of radiative hadronic background is fixed to 
the expectation N^adx = 1 -02 ± 0.14 for m^o < 4 GeV 
(m x < 2 GeV) where the IFR veto is applied, and 
to iV hadr = 6.6 ± 1.1 for fits in 4 < m A o < 8 GeV 
(2 < m x < 4 GeV) range. We do not observe a significant 
excess of events above the background, and set upper lim- 
its on B{Y{IS) -> jA a ) x B{A° -> invisible) (Fig. [2} and 
B(Y(1S) —> 7xx) (Fig. El). The limits are dominated by 
statistical uncertainties. The largest statistical fluctua- 
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FIG. 1: Projection plots from the fit with N s i s — onto (a,c,e) 
M lecoa and (b,d,f) M%. (a,b): low- mass region with IFR veto; 
(c,d): low-mass region without IFR veto; (e,f): high-mass re- 
gion. Overlaid is the fit with N s ig = (solid blue line), contin- 
uum background (black dashed line), radiative leptonic Y(1S) 
decays (green dash-dotted line), and (c,d) radiative hadronic 
Y(1S) decays or (e,f) r\ background (magenta dotted line). 




FIG. 2: 90% C.L. upper limits for B(Y(1S) 
invisible) . 
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FIG. 3: 90% C.L. upper limits for B(Y(1S) ->■ j X X)- 



tion, 2. Oct, is observed at m^o = 7.58 GeV [26j; we esti- 
mate the probability to see such a fluctuation anywhere 
in our dataset to be over 30%. 

In summary, we find no evidence for the single-photon 
decays T(1S) — > 7 + invisible, and set 90% C.L. up- 
per limits on B(Y(1S) ->■ jA°) x B(A° ->■ invisible) in 
the range (1.9-4.5) xl0~ 6 for < m A o < 8.0 GeV, (2.7- 
37)xl0~ 6 for 8 < m A o < 9.2 GeV, and scalar A . We 
limit B{Y{1S) -> jxx) in the range (0.5-24) xl0~ 5 at 
90% C.L. for < m x < 4.5 GeV, assuming the phase- 
space distribution of photons in this final state. Our 
results improve the existing limits by an order of magni- 
tude or more, and significantly constrain [26j light Higgs 
boson [13[ and light dark matter [8J models. 
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FIG. 4: Projection plots from the fit with m A o = 7.58 GeV (the most significant deviation from zero) to (a) Af roco n and 
(b) M x . Overlaid is the fit (solid blue line), signal contribution (solid red line), continuum background (black dashed line), 
radiative leptonic T(1S) decays (green dash-dotted line), and radiative hadronic T(1S) decays (magenta dotted line). The top 
plot show residuals in each bin, normalized by the bin error. The fit corresponds to B(T(1S) — > 7^4°) x B(A° — > invisible) = 
(3.2l^'g ± 1.0) x 10 -6 , where the first uncertainty is statistical and the second is systematic, and statistical significance of 2.0a. 
The probability to observe such a fluctuation anywhere in our dataset is over 30%. 
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FIG. 5: Upper limits on the product gr x \/B(A° — >• invisible) at 90% C.L. as a function of rn A o. The parameter g r is an 
effective coupling of the CP-odd Higgs A to bound state T(IS'); in NMSSM, gr — tan /3 cos 8 Fr , where cos8 is the fraction 
of non-singlet component in A , tan/3 is the ratio of Higgs vacuum expectation values, and Fr is the effective form-factor 
(including the QCD and QED corrections). The theoretically preferred region in NMSSM \\3{ is gr > 1- 



